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1. Introduction 

Radiative meson decays represent an important source of information on low- 
energy hadron physics, shedding light, for example, on the structure and properties 
of low- mass resonances, such as the /o(980). In particular, radiative decays to 
77 and 7/ can provide insights into the long standing problem of 77 — 17' mixing and 
probe the strange quark content of the light pseudoscalars [1, 2]. Radiative decays, 
not only raise interesting theoretical issues, but are an important focus of the data- 
taking by the KLOE experiment at the DA$NE 0-factory [3], where a large sample 
of decays will be collected, dramatically improving the experimental information 
already obtained by the VEPP— 2M groups at Novosibirsk [4]. 

This paper is devoted to the analysis of the radiative — > 777 and — > 7/7 
transitions. In contrast to the light vector meson case, where the u and are 
recognised as almost ideally mixed states with quark content of well defined flavour, 
77 — rj' mixing is still a much debated subject. The once conventional description was 
to adopt a single mixing angle in the octet-singlet flavour basis. Various attempts to 
estimate such an angle lead to results ranging from —10° to —20° [5]. More recently, 
Leutwyler et al. [6, 7] have shown that a consistent treatment of the 77 — rj' system 
requires the introduction of two mixing angles with a consequent redefinition of the 
particle decay constants. An equivalent description, as explained in more detail 
below, is obtained if the mixing basis is chosen to be the quark-flavour basis instead 
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of the octet-singlet one [8]. In such a scheme, it has been shown that a description 
in terms of a single mixing angle is quite reliable leading to predictions satisfying 
constraints from Chiral Perturbation Theory [6]. 

In this context a central role is played by the U(1)a anomaly. Since the flavour- 
singlet axial vector current is not conserved due to this anomaly, the 7/ meson cannot 
be identified as the ninth Goldstone boson. This crudely explains the fact that the 
7]' is much heavier than the other members of the pseudoscalar nonet. By combining 
chiral symmetry with the concepts of the large N c limit of QCD, Leutwyler[9] has 
extended Chiral Perturbation Theory from an expansion in the light quark masses 
and momenta to encompass powers of 1/N C . Then in the limit N c — > oo, the U(1)a 
anomaly vanishes and the rj' can formally be identified with the ninth Goldstone 
boson. In order to make the framework more predictive and closer to reality, cor- 
rection terms are added to the effective Lagrangian of the theory, expressing the 
deviation from the chiral limit for the decay constants and the masses of the light 
mesons. A Wess-Zumino- Witten (WZW) term describing the anomalous coupling to 
photons must also be added. Radiative decays provide additional information on 
the strength of this WZW term. 

In the following we analyse <fi radiative decays using QCD sum-rules [10] at lead- 
ing order in perturbative QCD. These sum-rules are widely recognised as a reliable 
technique for including the effects of non-perturbative QCD. In section 2 we survey 
possible 77 — 77' mixing schemes, with particular emphasis on the flavour basis mixing 
scheme developed by Feldmann et al. [8]. We give the relation between the parame- 
ters characterising such a scheme and those in the octet-singlet basis. In section 3, 
as a preliminary to our analysis of <p radiative decays, we compute the coupling of 
the strange pseudoscalar current both to the 77 and 7/ mesons using two-point QCD 
sum-rules. These provide key inputs into the three-point QCD sum-rule developed 
to estimate the decay widths T(<f> — > 7/7) and T(<f) — > 7/7), in section 4. The possible 
effect of 0(a s ) corrections is discussed at the end of this section. 

Though we often refer to the problem of mixing and we actually work with 
interpolating currents defined in the quark flavour basis, our results do not depend 
on any specific mixing scheme and so provide a genuinely mixing scheme independent 
set of QCD predictions. If then a particular flavour mixing scheme is adopted, our 
results can be translated into a prediction for one of the mixing angles. In section 
5 we again use two-point QCD sum-rules to compute the coupling of 77 and 77' to 
the strange and non-strange axial currents, identifying the results with the decay 
constants in the flavour basis mixing scheme. The results again allow us to estimate 
the mixing parameters in such a scheme. The results in sections 3 and 5 can be 
exploited to obtain an estimate of the contribution of the anomaly to the couplings 
computed in section 3. In section 6 we draw our conclusions. 
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2. On 77 — 77' Mixing 

Let us first recall the usual parametrization of 77 — rj mixing in the octet-singlet basis. 
We define current-particle matrix elements as 

<0\J^\P(p)>=t f P p, (i = 8,0; P = V,V'), (2.1) 

with Jf the SU(3)f octet axial vector current: 

Jl» = 7^=7 (m7 m 75M + d^ 5 d - 2s7 M 7 5 s) (2.2) 

and J° M the singlet current: 

J°n = 7^ (uj^u + d>y^ 5 d + 57^755) . (2.3) 

As already mentioned, two mixing angles, 6$ and 6q, are required [9] in order to treat 
mixing consistently. Accordingly, the couplings in (2.1) can be defined as follows: 

fl = fs cos # 8 /° = -f sin^o 

f* =f 8 sin9 8 /°=/ o cos0o- (2-4) 



Alternatively we can consider two independent axial vector currents with distinct 
quark flavour: 

J l» = si^s . (2.5) 



The couplings of the rj and rj' mesons to the currents (2.5) can be defined analogously 
to (2.1). The decay constants are written according to the following mixing pattern: 

f% = U COS 09 frj = ~fa Sin 0s 

frf = fq Sm 0<? fv' = fs C0S( Ps ■ (2.6) 

Though there are, of course, two angles in each basis, Feldmann [8] has shown 
that the mixing is specified quite accurately in terms of a single mixing angle, i.e. 
(fiq = (fis = 0, since \4> s — (p q \/ ((p s + (pq) 1, resulting in a much simpler framework. In 
this approximation, the states follow the same mixing pattern as the decay constants: 

\rj > — cos \i] q > — sin \rj s > 

1 7/ > = sin0 \r) q > +cos0 \i] s > (2.7) 

where \i] q > and \i] s > have a quark content defined by ideal mixing. We will refer 
to this simply as mixing in the quark-flavour basis, in order to distinguish it from 
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the previous one, which will be referred to as mixing in the octet-singlet basis. It 
is straightforward to obtain the relations between the parameters in the two mixing 
schemes: 

, a f q sin0 g - V2fs cos0 s f s sin0 s - ^2f q cos^ 

tan 9 8 = — p , tan 9 = p-^ 1 (2.8) 

f q cos <f) q + V2/ s sin S f s cos S + \/2f q sin g 

which become, for (p q = S = [8]: 

6 8 = (p — arctan ^"^y^^ #o = — arctan ^ "^y^" ^ ■ (^-9) 



Moreover: 



/ 8 2 = ^ 2 + ^ 2 + ^/J S sin(0 s -0 g ) 
/o 2 = + Ifs - ^fjs sin(0 s - 0,) . (2.10) 



In the following sections we derive the ingredients necessary for the description of the 
radiative — ^ 777 and — > 7/7 decays, without assuming any specific mixing frame- 
work. Then, in section 5, we will evaluate the couplings of the 77 and rj' to the axial 
vector current and estimate the parameters appearing in (2.6), obtaining a predic- 
tion for (pq, S . At the end, we shall comment on the accuracy of the approximation 



3. Two-point function for 77 and rf couplings to the pseu- 
doscalar current 

Let us consider the matrix element of the divergence of the axial- vector current: 

< 0|^J 5 > >=mtf. (3.1) 

As is well known, this divergence contains the axial-vector anomaly: 

d»Ji» = d^s^s) = 2m s si lb s + -^GG , (3.2) 

where G is the gluon field strength tensor and G its dual. This gives a relation 
between the matrix elements of the axial-vector current and of the pseudoscalar 
current: 



2m s < 0|sry 5 s|r7 (/) >= /Jo"iJ(0 - (° 
Let us call: 

< 0|sry 5 s|?7 >= A , (3.4) 



-^GG 
An 



V (>) ) ■ (3.3) 
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and compute this quantity by QCD sum-rules starting from the two-point correlator: 
T A (q 2 ) = i J d A xe^ x < 0\T[J s 5 (x) J 5 sf (0)]|0 > (3.5) 
where J| = si^^s. The correlator (3.5) is given by the dispersive representation: 

j roc q( g\ 

T A (q 2 ) = - ds + subtractions . (3.6) 

7T J Ami s - Q 

In the region of low values of s, the physical spectral density contains a 5— function 
term corresponding to the coupling of the rj to the pseudoscalar current. Picking 
up this contribution, we can write (dropping possible subtractions which we discuss 
later) : 

T ( 2\ A " 1 f°° A P kad ^ (1 7\ 

TA{q = — 2 + - / ds 2 ■ 3 - 7 

This corresponds to assuming that the contribution of higher resonances and contin- 
uum of states starts from an effective threshold s . On the other hand, the correlator 
Ta{q 2 ) can be computed in QCD by expanding the T-product in (3.5) by an Oper- 
ator Product Expansion (OPE) as the sum of a perturbative contribution plus non- 
perturbative terms which are proportional to vacuum expectation values of quark 
and gluon gauge-invariant operators of increasing dimension, the so called vacuum 
condensates. In practice, only a few condensates are included, the most important 
contributions coming from the dimension 3 < qq > and dimension 5 < qgaGq >. 
Here we follow such a prescription. 

In the QCD expression for the two-point correlator considered, the perturbative 
term can also be written dispersively, so that: 

t qcd,2^_ ds P \J. + d z <ss>+d b <sg(rGs >+..., (3.8) 

7T J Am 2 a s ~ Q 

where the spectral function p QCD and the coefficients ofa, d 5 can be computed in 
QCD. 

The next step consists in assuming quark-hadron duality, which amounts to 
assuming the physical and the perturbative spectral density are dual to each other, 
in the sense that they should give the same result when integrated appropriately 
above some s . This leads to the sum-rule: 

A 2 l r so p QCD (s) 

— - = -/ ds Y- + d 3 < ss > +d 5 < sgaGs > +... (3.9) 

m* — q z 7T J Ami s — Q 

This expression can be improved by applying to both sides of (3.9) a Borel transform, 
defined as follows: 
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where / is a generic function of Q 2 = —q 2 . The application of such a procedure to 
the sum-rules amounts to exploiting the following result: 



B 



-s/M 2 



[s + Q 



2\n 



(M 



2\n 



n - 1)1 



(3.11) 



where M 2 is known as the Borel parameter. This operation improves the conver- 
gence of the series in the OPE by factorials and, for suitably chosen values of M 2 , 
enhances the contribution of low lying states. Moreover, since the Borel transform 
of a polynomial vanishes, it is correct to neglect subtraction terms in (3.6), which 
are polynomials in q 2 . The final sum-rule reads: 



A 2 e m 



— m s e 
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m 2 






87T 2 J Am 2 
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< ss > ^ 
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4m 2 
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m; mi \ 1 ^ 1 1 mi 

W + W) + JP < ~ a!3aG! ' > V-2KP 



.(3.12) 
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In the numerical evaluation of (3.12) we use < ss >= 0.8 < qq >, < qq >= (— 0.24)' 
GeV 3 , < sgaGs >= 0.8 GeV 2 < ss >, m v = 0.548 GeV. 

The strange quark mass is chosen in 
the range m s = 0.125 - 0.140 GeV, ob- 
tained in the same QCD sum-rule frame- 
work [11]. The threshold is chosen be- 
low the T]' pole and varied between s = 
0.9 2 - 0.95 2 GeV 2 . 

Since the Borel parameter has no phys- 
ical meaning, we require that the result 
does not depend on it. This is achieved 
by finding a "stability window", i.e. an 
interval of values of M 2 , where the out- 
come of the sum-rule is almost indepen- 
dent on M 2 . Such a window is usually 
sought in a restricted interval of values of 
the Borel parameter chosen by requiring 
that the perturbative contribution is at 

least 20 % of the continuum (which corresponds to considering the integral in the 
perturbative term up to infinity rather than up to Sq), which produces an upper 
bound on the Borel parameter: M 2 < 1 GeV 2 . Additionally requiring that the per- 
turbative term is greater than the non-perturbative contribution, the lower bound: 
M 2 > 0.5 GeV 2 , is obtained. Then, the stability window for M 2 in [0.8, 1] GeV 2 can 
be selected. 



M 2 (GeV 2 ) 



Figure 1: Coupling of the r/ to the pseu- 
doscalar current as a function of the Borel 
parameter M, for m s = 0.133 GeV. The 
solid curve corresponds to the higher thresh- 
old so = (0.95 GeV) 2 , the dashed curve cor- 
responds to sq = (0.9 GeV) 2 . 
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In figure 1 we plot the sum-rule (3.12) for m s = 0.133 GeV, which corresponds to 
the central value of the range of values adopted in the analysis. Taking into account 
the uncertainty on m s , we obtain: 

\A\ = (0.115 ± 0.004) GeV 2 . (3.13) 

Some comments are in order on the accuracy of the result (3.13). This has been 
obtained at leading order in perturbative QCD, as with all the results presented in 
this paper. Consequently, the uncertainty affecting the determination (3.13) should 
be taken modulo the neglect of a s corrections, the role of which we comment on later. 

Another source of uncertainty is linked to the choice of the strange quark mass. 
It should be stressed that the value of this parameter is quite controversial. On the 
one hand, lattice determinations seem to point towards lower values of m s (for recent 
reviews see e.g. [12]); on the other hand, results obtained in other approaches indicate 
higher values [13]. As for QCD sum rules, the result in [11] exploits an accurate 
determination of the hadronic spectral function based on experimental information 
on the Kn system including a non-resonant component in addition to the resonances 
in the 1 = 1/2 channel. In [11] it has been shown that the effect of this non- 
resonant contribution is a reduction in the spectral function, with a consequent 
lowering of the value of m s with respect to previous QCD sum rule determinations 
[14]. We have therefore chosen to adopt the result of [11] at first for consistency, 
i.e. using a result obtained by the same technique, but also because such a value for 
m s falls in the middle of the existing range. Moreover, the value fh s (l GeV)=0.125 
GeV obtained in [11] could be considered as a lower bound on this parameter, since 
further experimental information could be added to improve the sum rule further. 
Consistently, we have used the range for rh s of 0.125-0.140 GeV quoted above. The 
result (3.13) turns out to be quite stable. Indeed, we have explicitly checked that 
using still higher values for m s (up to 0.160 GeV) would produce little change. 

Let us now consider: 

< 0\si-f 5 s\r]' >= A' . (3.14) 

An analogous calculation gives: 



(A ) e + A 2 e TP = — / ds s \ 1 -e 



<ss> 1-t^| + t^t +t7^< sgaGs > 1 1 



M 2 M 4 / M 2 * \ 2.1/ 2 



(3.15) 



where we have raised the effective threshold up to s' in such a way as to pick up the 
i]' pole too: Sq = (1.44, 1.55) GeV 2 . Using = 0.958 GeV and fixing the stability 
window for M 2 to be [1.2,2] GeV 2 , we obtain: 

\A'\ = (0.151 ± 0.015) GeV 2 . (3.16) 
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We shall use the results (3.13), (3.16) in the next section. Though we cannot actually 
establish the sign of A, A' from the sum rule, we assume that A ■ A' > 0. 



4. Radiative — > 777 and — > 7/7 decays 

Having found the key matrix elements A and A' of (3.13), (3.16), we now consider 
the three-point functions defined by 

< 77(^)1*7^10(21,61) >= F(q 2 ) e^ s (q 1 Uq 2 ) p (e 1 ) s (4.1) 

(? = Qi ~ 92)- In order to compute the 777 decay, we need the coupling (7 = 
— |F(0), obtained for a real photon coupling to a strange quark. We consider the 
three-point function: 

n^(g 2 ,g 2 2 ,g 2 ) = i 2 J d A x d A y < 0\T[J°(y)J u (0)J,(x)]\0 > (4.2) 

where J| has been defined above and J v = 57^ s is the vector current. The correlator 
(4.2) can be written as: 

IW^ 2 , ql q 2 ) = Tl(ql q 2 2 , q 2 ) e^( qi ) a (q 2 f (4.3) 

and a QCD sum-rule can be built up for the structure U(q 2 ,ql,q 2 ). The method 
closely follows the one described for the two-point sum-rule. We assume H.(q 2 , q%, q 2 ) 
obeys a dispersion relation in both the variables ql,q 2 : 

nte 2 , ql q 2 ) = ~J dsi / ds 2 P^ f^ , (4.4) 
tt 2 J J (si - qi){s 2 - q 2 ) 

with possible subtractions. Such a representation is true at each order in perturbation 
theory and, as is standard in QCD sum rule analyses, it is assumed to hold in 
general. In this case the spectral function contains, for low values of si, s 2 , a double 
5— function corresponding to the transition — > 77. Extracting this contribution, we 
can write: 

rrf 2 2 2^ AF{q 2 )m 4> f 4> 1 /•<*> /■«> p fead (si, gg, g 2 ) ,x 

n(9i,9 2 >9 ) = 7— 2 2vHr 2T + ^/ rfs i / ^27 2Y7 2T ' 45 

(mj-qi){ml-q$) tt 2 J S01 Js 02 ( Sl - qf)(s 2 - ql) 

where subtractions are neglected as later they will vanish on taking a Borel transform. 
The parameter A appearing in the previous equation is just the coupling of the 77 
to the pseudoscalar current, computed in section 3. Deriving an OPE-based QCD 
expansion for II for large and negative q 2 , q\ and q 2 , one can write: 

1 r°° . r°° . p QCD (si,so,a 2 ) 



tt, 2 2 2n f°° j f°° j P^~{Sl,S 2 ,q 

n(gi,g 2 >9 ) = — / d- s i / ,« s 2 

IX J Ami J Ami 



si - q 2 ){s2 - ql) 
+ c 3 < ss > +c 5 < sgaGs > +... . (4.6) 
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Invoking quark-hadron global duality as before, we arrive at the sum-rule: 

[ dsids 2 - 

7T JD 



AF(q 2 )m^U = l_r p QCD ( Sl , s 2 , g 2 ) 

(ml - q 2 ){ml - q\) ix 2 Jd 12 (si - qf)(s 2 - q\) 



+ c 3 < ss > +c 5 < sgaGs > +... . (4.7) 

where the domain D should now also satisfy the kinematical constraints specified 
below. After a double Borel transform in the variables —q\ and — q%, we obtain: 



AF(q 2 )m^ = e^e^\ \ ds x / ds 2 e ^ e ^ 8 

[ J J 7i 2 ^X(s 1 ,s 2 ,q 2 ) 



+e M i e M 2 



/ m 2 m 2 mi m\ m 2 (2m 2 — g 2 V 
<ss> 2 - — ^ - — ^ + — s -r + — s -r + ' 



Ml Ml Mf Mi M 2 M. 



2 



(I 2 m 2 m 2 (2g 2 -3m 2 ) 



6M 2 ' 3M| 2Mf 2M 2 4 ' ZM%M$ Jlj^ 

The integration domain D over the variables s±, s 2 depends on the value of q 2 and is 
given by D = D 1 U D 2 where: 

• (~q 2 ) > $02 - 4m 2 

Di- fa)- < s 2 < s 02 4m 2 < si < s i 

• i-q 2 ) < SQ2 - 4m 2 

D 2 ■ (s 2 )- < s 2 < (s 2 )+ 4m 2 < si < (si)_ 

(s 2 )- < S2 < S02 < si < s i (4.9) 

with: 



2m 2 <? 2 + (2m 2 - o 2 )si ± \ Siq 2 (q 2 - 4m 2 ) (si - Am 2 ) 
(s 2 )± = — ^_J^ ^ =i (4.10) 



s 



2m 2 s q 2 + (2m 2 -q 2 )s 02 ±Js 02 q 2 (q 2 -Am 2 )(s 02 -Am 2 ) 
(«i) ± = ^ • (4.H) 

Since we consider the form-factor F(q 2 ) for arbitrary negative values of q 2 , we 
could perform a double Borel transform in the two variables Q\ = —q\ and Q\ = —q 2 , 
which allows us to remove single poles in the si and s 2 channels ( "parasitic" terms) 
from the sum-rule. Our procedure is therefore to compute the form-factor F(q 2 ) 
and then to extrapolate the result to q 2 = 0. Strictly speaking, since we only know 
the magnitude of A in (4.7), it is the modulus of F(q 2 ) that is determined. In the 
numerical analysis we use: m^ = 1.02 GeV, = 0.234 GeV (obtained from the 
experimental datum on the decay to e + e _ [15]). We compute the result for two 
values of the <fi threshold: Sqi = 1.8, 1.9 GeV 2 . Sq 2 coincides with the f] threshold 
chosen as we did for the two point function in section 3. 
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The outcome of the sum rule is de- 
picted in figure 2. Varying all the param- 
eters entering in the sum rule we obtain 
a region delimited by the dashed and the 
dotted curves in this figure: they corre- 
spond to the set of parameters giving the 
highest and the lowest curve analytically 
defined by (4.8). The resulting form fac- 
tor shows a behaviour in — q 2 in the region 
(-0.8, -0.2) GeV 2 which can be fitted by 
a parabolic function. The extrapolation 
to q 2 = gives: 

\ g \ = EM. = (0.66 ± 0.06) GeV- 1 . 

(4.12) 

The central value in (4.12) corresponds 
to the extrapolation of the solid line in 
figure 2. The uncertainty range in \g\, as 
obtained by the extrapolation procedure, 
is displayed in the same figure. 

We can now use this result to com- 
pute the — > 777 decay width: 




Figure 2: Form factor F{q 2 ) obtained 
varying the input parameters in the sum 
rule (4.8). The isolated point on the right is 
the result of an extrapolation. The extrap- 
olation of the solid curve gives the central 
point on the right, corresponding to the re- 
sult (4.12). 



T(0 777) 



24 



m 



m~ 



777,4, 



and, using T(0) = 4.43 MeV [15], we obtain: 



B((f> -> 777) = (1.15 ±0.2) 



777) 



(4.13) 



(4.14) 



1.18 ± 



which compares favourably with the experimental outcome: B((p 
0.03 ±0.06)% [4]. 

We can extend the previous analysis to the channel with the 77' in the final state. 
The derivation of the sum-rule is straightforward: 



A' F V W/*e M 2 + A F{q 2 ) m(t> Ue 

3m s 




dsi / ds 2 e M i e 



s 2 



7T 2 ^X(s 1 ,s 2 ,q 2 ) 



+ e 



'1 e 



<ss> 2 



mz 



M 2 



2 4 4 

mz mz ml 
— - H H + 

M 2 Mf Mi 



m 2 {2m 2 — q 2 
Ml Ml 
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+ < ?9<Tfe > (wF + 5M? " 5M? " 2Ag + 3«?Ag J J} (415) 

The integration region is the same as before with the substitution: s 02 — > s' 02 = s' 
(as in the two point sum-rule of section 3). We then obtain: 

\ g '\ = E!M = (i.o ± 0.2) GeV- 1 (4.16) 

o 

which yields 

B((j)^ rfi) = (1.18 ±0.4) 1(T 4 . (4.17) 

The experimental datum is: B{<j> -> 7/7) = (0.82lg;S ± 0.11) 10~ 4 [4], completely 
compatible with our result (4.17). 

These results have been derived without including QCD radiative corrections. In 
principle, each term in the OPE could be computed as an expansion in powers of a s , 
supplementing the non-perturbative expansion with short-distance corrections. This 
would display the correct scale and scheme dependence for the hadronic quantities, 
such as the coupling of the 77 and 7/ to the currents considered in our analysis. The 
calculation of QCD corrections is a difficult task well beyond the scope of the present 
paper. However, we would like to comment on the possible role of such terms. 

As far as the two point sum rule is concerned, 0(a s ) contributions have been 
computed [24]. At a typical scale /i — 1 GeV, these corrections are sizeable, indicating 
that still higher orders may also be important. On the other hand, the main goal 
of the present analysis is the computation of radiative decays. These results are 
obtained from the ratio of three point to two point sum rules. The most reliable 
procedure in this case is to compute consistently the three point and two point 
correlators at the same order in a s . The uncertainty due to the neglect of higher 
order corrections should be reduced due to a cancellation in the ratio. 

This expectation is fulfilled, for example, in the calculation using QCD sum rules 
of the Isgur-Wise function [25], describing in the heavy quark limit the B — > 
semileptonic transitions. In this case, though the 0(a s ) corrections are large for the 
two point sum rules [26], explicit calculation of the three point correlator shows the 
expected cancellation, i.e. the modest role of radiative corrections in the outcome. 
In this particular case, the result is expected, at least at the zero recoil point, by 
symmetry requirements. However, this cancellation works in more general situations, 
such as the one presented in [27], where the universal form factor describing the 
B transitions to orbitally excited charmed mesons was computed at order a s by 
the same method. Again, although the two point correlator received important 
corrections, the ratio ot three to two point functions is quite stable. 

In the light of this discussion, our results for the decay constants, eqs. (3. 13), (3. 16), 
should be considered to be estimates, the uncertainties in which do not take into ac- 
count possibly sizeable radiative corrections. On the other hand, the outcome for 
radiative decays, eqs. (4.14), (4.17), should be viewed as much more accurate. 
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Indeed, the predictions for the branching ratios in (4.14) and (4.17) are the major 
results of this paper. They are quite independent of any mixing scheme for the 77 
and rj'. Nevertheless, adopting the mixing scheme in the flavour basis described in 
section 2, it is possible to derive the relation: 

« = f^ = (4^4T^. (4.i8) 

from which we get <p s = (34±g)°. The experimental ratio would give: <p s = (39.O1L5 5) . 

As mentioned in the introduction, the results obtained can, in principle, provide 
us with information about the magnitude of the WZW term, which represents an 
OZI-rule violating contribution to the effective lagrangian. The strength of this term 
is parametrized by a constant A 3 and is determined by the values of the couplings g 
and g' . For example, in [5] it is found: 

_ 3m / 1 cos 4> s sin <p v l sin0 g A 3 sin6> 8 \ , > 

where (fiv is the mixing angle in the 0— u system. We assume this formula to estimate 
the size of A 3 . Unfortunately, the combined effect of the uncertainties affecting the 
parameters entering (4.19) allows us no more definite conclusion than A 3 ~ 2 ± 4. 
Precise measurements of <j) radiative decays to 77 and 7/ at KLOE [3] will do better. 

Let us now compare our results with previous determinations. In ref. [16] the 
chiral anomaly prediction at q 2 = and the vector meson dominance are exploited 
to derive the couplings g = g^ and g' = g^'-y. A single angle mixing framework 
in the octet-singlet basis is assumed and the corresponding coupling constants fo, f$ 
are derived from the experimental data on the decays 77 — > 77 and 7/ — > 77 and used 
as an input to derive g and g' as a function of the mixing angle. In [8] this approach 
is extended to the quark-flavour mixing scheme with the result: g = 0.78 GeV -1 and 
g' = 0.95 GeV -1 . On the other hand, an energy-dependent mixing scheme is adopted 
in [17], with the result: g = (0.73 ± 0.06) GeV -1 and g' = (0.83 ± 0.06) GeV -1 . 
Alternatively, Ref. [18] exploits the hidden local symmetry approach, together with 
the inclusion of various SU(3) symmetry breaking [19, 20] to obtain g — g' — 0.70 
GeV" 1 . As can be observed, the various approaches seem to agree quite well for the 
decay with the 77 in the final state, while the results have a larger spread in the case 
of the rj'. For a more comprehensive survey of results we refer to [5]. 



5. 77 and rf couplings to axial- vector currents using two-point 
sum-rules 

Our study of <fi radiative decays to 77 and 7/ require no knowledge of 77 — rj' mixing, 
only the coupling to strange quarks is needed. However, using similar techniques, 
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we can investigate their decay constants in both the strange and non-strange sectors 
and so deduce the mixing pattern within errors. This is the purpose of this section. 
We begin by considering the correlator 



(5.1) 



IW) =i d*x e«« < 0|T[^(x) J^(0)]|0 > 



Following the procedure already outlined above, we obtain the sum-rule: 



eM 1 



1 f s ° , » vert , , 2m s 



< ss > e m 1 



where: 



■ert i 



1 L 4m 2 s 
An V s 



2m 2 + s 



(5.2) 



(5.3) 



(in this case the d = 5 contribution vanishes). 

The allowed range for the Borel pa- 
rameter, obtained according to the above 
criteria, is: 0.35 GeV 2 < M 2 < 3.5 GeV 2 , 9 
and the further stability window is found $ 
in [2, 3.5] GeV 2 . The result is depicted in 
figure 3, for the value m s = 0.133 GeV. 
Taking into account the uncertainty in 
m s too we get: 




M (GeV 



f° = (0.13 ± 0.01) GeV . 



(5.4) 



Figure 3: Constant /* ELS £1 function of the 
Borel parameter M for m s = 0.133 GeV. 
The solid curve corresponds to so=(0.95 
GeV) 2 , the dashed one to s =(0.9 GeV) 2 . 



In order to determine /*/ one has to 
repeat the previous calculation more or 
less exactly, raising the threshold above 
the rf mass and considering the pole contribution of the i] on the hadronic side of 
the sum-rule. The result is: 



+ (/, 



t\2 * 



2m* 



71 J Ami 



dse 'm 1 ff (s) + —j < ss > e m 



(5.5) 



where p peri (s) is the same as before. In the numerical analysis, s' is varied in the 
range (1.2 - 1.25 GeV) 2 . The selected stability window for M 2 is [2,5] GeV 2 . We 
obtain: 

f s v , = (0.12 ± 0.02) GeV . (5.6) 

If we now use the current J| M in the correlator, instead of J| , and we set the up 
and down quark masses to zero, we can obtain The stability window is found for 
M 2 in the range [2, 4] GeV 2 with the result 



m 

J V 



(0.144 ±0.004) GeV . 



(5.7) 
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Raising the threshold, we can also evaluate with the result: 

= (0.125 ± 0.015) GeV . (5.8) 

As with the results obtained in section 3 from two-point sum rules, these also repre- 
sent estimates, derived without radiative QCD corrections. 

We can now use the set of results (5.4), (5.6), (5.7), (5.8) to estimate all the 
mixing parameters appearing in (2.6). From the relation f^/f^ = tan0 s , one has 
<p s = (46.6° ± 7°). Since /* = f s sm<p s \ a prediction for f s follows: f s = (0.178 ± 
0.004) GeV, the central value corresponding to f s = 1.345/^, with f w = 0.132 GeV. 

Using the relation: fi/ = tan0 9 , we find 4> q = (41° ± 4°). We can now derive 
a prediction for f q , using = f q cos0 g . We obtain f q = (0.19 ± 0.015) GeV, the 
central value of which corresponds to f q ~ lA4f n . Let us observe that our results 
correspond to \<p s — 4> q \/(4> s + <p q ) — 0.065, which confirms the relation put forward 
in [8] that this ratio should be much less than 1. 

If we now turn to the scheme with two mixing angles in the octet-singlet basis, 
we could exploit the relations (2.8)-(2.10) to obtain: 

# 8 ~ -8.4° # - -13.8° 

/ 8 ~ 1.44 U fo - 1-35 f n . (5.9) 

Previous determinations of the parameters calculated above range over large 
intervals, expecially for those corresponding to (5.9), i.e. in the octet-singlet mixing 
scheme. For a comprehensive collection of previous results we again refer to [5]. We 
only observe that our results for <f) q , <p s are in pretty good agreement with those in refs. 
[7, 18, 21]. Our outcome for f s also agrees quite well with most of previous results 
[5], while the result for f q seems somewhat larger than previous determinations. 

We can now exploit the results obtained in this section, i.e. the values in (5.4), 
(5.6), together with the predictions (3.13) and (3.16), to derive the following matrix 
elements from (3.3): 

rj) = (0.008 ± 0.004) GeV 3 , 

rj') = (0.072 ± 0.025) GeV 3 . (5.10) 

Both values in (5.10) are close to the naive quark model calculation of Novikov et 
al. [22], particularly that for the r/, which is f n m^/y/S = 0.070 GeV 3 . Our result 
for the rj is somewhat smaller than the one found in [22]: f n = 0.016 GeV 3 . 

However, their simple quark model result does not take into account SU (3) p breaking 
corrections. In contrast, ours does. Thus, the matrix elements of (5.10) are important 
for the investigation of the structure of the r\ and r/ and their possible glue content 
[22, 23]. 

1 This relation is to be considered in terms of absolute values, since the sum-rule gives access 
only to {f s q ) 2 1 and therefore does not allow the sign of /* to be determined. 
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6. Conclusions 

We have analysed radiative — > 777, <fi — > 7/7 decays using QCD sum-rules. This 
analysis required a preliminary calculation of the couplings of the pseudoscalar cur- 
rent to the r) and 7/. The sum-rules are derived without any assumption about f] — i] f 
mixing. Though we use only lowest order in QCD perturbation theory, potentially 
large higher order corrections are expected to cancel between the three and two point 
correlators, to give results that are in good agreement with the available experimen- 
tal data on the r\ channel, and are compatible with the Novosibirsk datum for the 
7/ case. Since the uncertainty in the latter is large, the last word is still left to the 
experimental improvement at DA$NE, for instance. 

We have also discussed the issue of 77 — 7/ mixing, giving predictions for the 
parameters describing such mixing in a quark-flavour basis scheme. We observe that 
the two angles required in such a scheme are quite close to each other. The existing 
spread of results gives us confidence that new experimental information will shed 
light on this sector of low-energy physics too. 
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